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ABSTRACT 



Optimal control theory is used to develop three automatic 
control systems for the longitudinal motion of a semisubmerged 
ship. A linearized mathematical model of the ship motion is 
used and the control problem is treated as a linear regulator. 
Simulations of the ship's longitudinal motions, utilizing the 
three control systems, are compared for various sea conditions. 
It is concluded that if the wave forces and moments are known 
or are estimated as functions of time, a suboptimal controller 
is the best controller. However, if a frequency domain approach 
is used to estimate the wave forces and moments, the complexity 
of a suboptimal controller approaches that of an optimal con- 
troller and the difference between the two is negligible. 
Preliminary work on three techniques for estimating the forces 
and moments is presented. 
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I. 



INTRODUCTION 



In 1968 the need arose for a small, relatively inexpensive 
ship which could travel as fast as large Naval vessels and 
provide support for one of the Naval Undersea Research and 
Development Center (NUC) projects — an unmanned undersea 
vehicle. The ensuing investigation conducted by NUC resulted 
in the semisubmerged ship concept shown in Fig. 1 [Ref. 1] . 

Referring to Fig. 1, the concept is described as follows: 

a. The two streamlined hulls provide the primary 
buoyancy and have less drag and experience smaller wave- 
induced forces than conventional ship hulls. 

b. A propeller is placed at the end of each hull 
where efficiency is increased by utilizing the boundary layer 
inflow. 

c. Two vertical, surface-piercing bouyant struts are 
attached to each hull, providing hydrostatic stability in pitch 
and heave due to the waterplane areas and spacing. Hydro- 
dynamic stability in yaw is obtained from the hydrofoil behavior 
of the struts. 

d. The rudders are placed in the propeller slipstreams 
which provides additional rudder force and allows the rudders 
to induce forces on the aft struts. 

e. A horizontal tail surface provides pitch damping as 
well as pitch stability at high speeds. 

f. The superstructure, aside from providing the lormal 
equipment and personnel spaces, provides the juncture fo_ the 
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Figure 1. Semisubraerged Ship Concept. 



rigid attachment of the two hull-strut systems. This rigid 
system provides stability in roll. 

g. The elevator control surfaces incorporated in the 
horizontal tail surface allow a means of controlling the ship's 
response to waves when used in conjunction with 

h. the canard control surfaces which are mounted on 
the forward part of the submerged hulls. 

The purpose of this investigation was to consider various 
ways of controlling the elevator and canard surfaces and make 
recommendations concerning the method of control and areas 
requiring further investigation. 



9 



II. PROBLEM FORMULATION 



A. NON-DIMENSIONALITY 

The response of a ship to forces — environmental and/or 
self-induced by control surface deflections — can be deter- 
mined from a system of equations which describe the ship. If 
the particular ship's displacement is changed, but the hull 
configuration is retained, it is necessary to determine a new 
system of equations in order to determine the ship's response. 

A non-dimensional system of equations however, will allow the 
response to be determined from this one system of equations 
and scaled to any particular displacement. The system of equa- 
tions in this' investigation are written in a non-dimensional 
form based on Froude number scaling. 

Froude number = ■ — (2.1) 

/g V 1/3 

where 

U = ship's forward velocity 
g = acceleration of gravity 
V = ship's displacement volume 

The dimensional values can be determined by multiplying 
the non-dimensional values by their appropriate factors as 
given in Table I. The symbol p indicates water ' density . The 
characteristic length, l c , is the cube root of the displace- 
ment volume, V. The ship's displacement weight is denoted by 
the symbol A. 
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TABLE I 



CONVERSION FACTORS 
(Non-dimensional to Dimensional) 



Variable 


Multiplication 

Factor 


Length 




1 = v 1 / 3 

c 


Area 




_ v 3 / 3 

• c 


Volume 




O U> 
1! 

<3 


Linear Velocity 




/g 1 
^ c 


Time 




/T c /g 


Angular Rate, Frequency 


' 571 ; 


Linear Acceleration 


g 


Angular Acceleration 


9^0 


Force 




pgV = A 



B. SEMISUBMERGED SHIP FORM USED FOR THEORETICAL ANALYSIS 

The geometric form used in this discussion and in investiga- 
tions conducted by Naval Undersea Center (NUC) , San Diego, 
California, is shown in Fig. 2. The dimensions are given in 

terms of the characteristic length, 1 . Figure 3 will aid in 

c 

interpreting the results for a specific ship size. 

C. COORDINATE SYSTEM 

A body moving in a fluid can move in all six degrees of 
freedom of motion — translation along each of three orthogonal 
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Figure 2. Semi submerged Ship Form Used for Theoretical 
Analysis [Ref. 3] . 
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Displacement > A , long tons 
(Ship length ~ 5 V /} ) 



Figure 3. Basic Scaling Parameters as a Function of Ship 
Displacement [Ref. 3]. 
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axes and rotation about each of the three axes. The choice 
of a coordinate system should result in one in which the 
equations of motion can be most conveniently analyzed. 

The ship under investigation has a plane of symmetry (as have 
most ships in existence), which is defined by the ship's 
centerline and a line perpendicular to the ship's deck. Two 
of the axes are chosen to lie in this plane of symmetry with 
the remaining axis being orthogonal to the plane. This choice 
simplifies the expressions for the hydrodynamic forces through 
symmetry and simplifies the equations of motion because the 
axes are parallel to the principal axes of inertia [Ref. 2] . 

If the axes are located such that their origin coincides with 
the ship's center of mass, the axes become the principal axes 
of inertia. 

The ship's position is given in an Earth-fixed right-hand 
cartesian coordinate system. The attitude of the ship can be 
given by an Euler angle transformation -from axes parallel to 
the Earth-fixed axis system to the ship-fixed axis system as 
chosen above. Figure 4 shows the relationship between the two 
systems . 

1. X - Axis 
s 

The X g - axis is the ship-fixed longitudinal axis in 
the plane of symmetry with the positive direction being toward 
the bow. It is parallel to the ship's keel. 
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Figure 4. Coordinate System 

2 . Y - Axis 

s 

The Y g •- axis is the ship-fixed transverse axis per- 
pendicular to the plane of symmetry. The positive direction 
is measured to starboard. 

3. Z - Axis 

s 

The Z - axis is the "vertical axis" in the plane of 
s r 

symmetry. The positive direction is measured downward toward 
the keel. 

4 . Angles of Rotation 

The axis system that has been chosen forms a consistent 
right-handed coordinate system. If <(> is the roll angle, 9 is 
the pitch angle and is the yaw angle, then the positive rota- 
tions are as indicated in Fig. 5. 

D. LONGITUDINAL RESPONSE PROBLEM 

Previous investigations of the semisubmerged ship ress nse 
[Ref. 3) have indicated that the pitch and heave motions in 
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Figure 5. Angles of Rotation 

regular head waves are quite small. However, the ship's 
natural longitudinal response to following waves having a 
height- to-wavelength ratio of 0.05 has been found to be greater 
than desired for wavelengths greater than one ship length. In 
particular, with a non-dimensional speed of 1.65 (high speed), 
the heave response to a wavelength of 10 (height of 0.5) has 
been found to be greater than 1.25. The values for a 2000-ton- 
displacement ship (refer to Fig. 2 for dimensions) would be 

wavelength, X = 410 ft. 

wave height, h = 20.5 ft. 

heave , z = 51.25 ft. 

This' indicates that the ship is thrown out of the water, then 
slammed back into the water with the upper platform submerging. 
This is an intolerable situation and calls for the design ~>£ a 
control system to reduce the ship's natural response. 
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III. SIMULATION MODEL 



A. LINEARIZED EQUATIONS OF MOTION 

The force and moment equations for the ship are written 
in the ship-fixed axis system previously described. If it is 
assumed that the ship is running at a constant forward veloc- 
ity, U, and that all other motion variables (velocities, angular 
rates and attitude angles) make only small perturbations about 
straight equilibrium running, then the equations of motion can 
be linearized and separated into two groups: (1) the lateral- 

direction equations, and (2) the longitudinal-direction equa- 
tions (Ref. 2] . 

The longitudinal-direction equations of motion [Ref. 3] 
are written below in non-dimensional form. The vertical force 
is given by Eq. (3.1), the pitching moment by Eq. (3.2) and the 
kinematic relationships by Eqs. (3.3). The notation is de- 
fined in Table II. 



w'-Uq 



A A 



1/2 U [C w+C q 
w q 

+C 6 +C 6 ] - A tT z 
Zg 6 z ^ c w L 

°e °c 

+a T x e+z- ,w' 

wL wL w ' 



A A A 

+ Z- ,q' + Z 
q ^ w 



(3.1) 
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TABLE II 



DEFINITION OF NOTATION FOR EQUATIONS OF MOTION 




6 

/\ 

Z 




□ 

T 





/S 




M/ (1/2 p U 2 V) 



Non-dimensional wave forces along the X g axis 
Non-dimensional wave pitching moment about the Y g axis 
Non-dimensional water plane area 

Non-dimensional X coordinate of the centroid of water 

s 

plane area 

Non-dimensional metacentric height in pitch 
Non-dimensional radius of gyration about the pitch axis 
Y s 

Pitch angle 

Non-dimensional position of ship center of mass in 
Earth-fixed coordinates 

Elevator and canard surface deflection angles 

Velocity component along X and Z axes respectively 

b b 

Angular rate about Y g axis 

Force along Z g axis 

Moment about Y axis 
s 

Non-dimensional form of □ 

1/3 1/2 

Non-dimensional time, t(g/V ' ) ' 

da/dx 



M 



3 9 3 

3F L ' "3a 

Z/ (1/2 pU 2 V 2/3 



C 

3d z 
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2~ 

p q' 

y M 



= 1/2 U [C w+c q 

m^ 

w q 



+C <5e+C 6 ] - z 0 

m ■" mp 



m 



m 6 c ' 
c 



(3.2) 



A /V ^ ^ 

+A T X T z + M~ . w ' 
wL wL w ' 

+M^ , q ' + M 
q ' ^ w 



0' = q 

A A A 

z ' = — U 0 + w 



(3.3) 



The equations of motion can be rewritten in a form where 

A 

they are functions of pitch (0) , heave (z) , their respective 
first and second derivatives, control surface deflections and 
wave forces and moments by utilizing Eqs . (3.3). These equa- 

tions are shown below with the coefficients replaced by their 
numerical values as determined by previous investigations 
[Ref. 3). 

z" = 1/2 U [-4 . 303 ( z 1 +U0) 

-4 . 7630 1 -1 . 50U 6 

! (3.4) 

-0.57U 6 ] -0 . 60 8z 

-0.06470-1.075 (z " +U0 ' ) 

-0.3640" + Z 

w 

1.4120" = 1/2 U [-4 .763 (z+U0) 

-15.080 ’-2. 85U 6 

© 

+0.90U 6 )-l.O450 (3,5) 

c 

-0.0647z-0.364(z" +U0' ) 

-2.1420" +M 

w 
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B. TRANSFER FUNCTIONS 

Reference [3] gives the response of the ship in pitch and 
heave resulting from 0.1-radian control surface deflections. 

These responses can be determined from the transfer functions 
relating the responses to the control surface deflections. In 
order to determine if the model used in this investigation was 
equivalent to the one used by NUC, San Diego, the transfer 
functions were derived as shown below and the responses to the 
control surface deflections were obtained. 

The linear Eqs. (3.4) and (3.5) can be Laplace transformed; 
the result of doing this is shown below with the assumption made 
that the initial conditions are zero. The numerical values of 
the coefficients a^ and b^, i=l,...,8, are given in Table III. 

A 2 A A 2 

z(s)s = a^z (s) s+a 2 z (s) +a^6 (s) s 

+a^0 (s) s+a^S (s) +a g<$ e (3.6) 

+a 7 S c +a 8 S w 

0(s)s 2 = b^0 (s) s+b 2 0 (s) tb^z (s) s 2 

+b. z (s) s+b c z (s) +b,6 (3.7) 

4 o be 

+b_6 +b 0 M 
7 c 8 w 

One method of determining transfer functions is to draw 
the signal flow graph corresponding to Eqs. (3.6) and (3.7) and 
then apply Mason's Gain Rule to the resulting graph [Ref. 4]. 
Referring to Fig. 6 and Table III for the appropriate coefficient 
values, it can be seen that the transfer functions will we the 
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TABLE III 



Coefficient 


Value 


Coefficient 


Value 


a l 


-1.0369 U 


b l 


-2.2240 U 


a 2 


-0.2930 


b 2 


-0.6701 U 2 -0 .2940 


a 3 


-0.1754 


b 3 


-0.1024 


a 4 


-1.6658 U 


b 4 


-0.6701 U 




'-2 






a 5 


-1.0369 U -0.0312 


b 5 


-0.0182 


a 6 


-0.3614 U 2 


b 6 


-0.4010 U 2 




~ 2 




^ 2 


a 7 


-0.1373 U 


b 7 


+0.1266 U 


a 8 


+0.4819 


b 8 


+0.2814 



A 




A 



Z($) 



0CS) 



Figure 6. Signal Flow Graph. 
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ship's non-dimensional speed, U, as the determining parameter. 
The resulting transfer functions are of the form 



X(s) = 

YlST 



n 2 s +n 1 s+n q 



D 4 s 4 +D 3 s 3 +D 2 s 2 +D 1 s+D Q 



(3.8) 



where the coefficients N^, i=0,l,2, and , j=0,l,...,4 are 
given in Table IV. 

Reference [3] investigates the ship's response for two 
distinct speeds, u=l . 0 and 1.65. Referring to Fig. 3, it can 
be seen that for a 2000-ton-displacement ship, these values 
correspond to speeds of 21.8 knots and 36 knots respectively. 

The transfer functions for these two speeds are given in Table 
V, and the responses resulting from 0.1-radian control surface 
deflections were found to agree with those contained in Ref. [3] 



C . COMPUTER MODEL 

The time histories of the ship's responses, pitch, pitch 
velocity, heave and heave velocity, can be found by using a 
computer. A choice, however, must be made whether to use analog 
or digital simulation for the ship model. Both models are given 
below. 

1 . Analog Model 

An appropriate model for analog simulation can be easily 
achieved by normalizing Eq. (3.8) with repect to D^. The result 
ing model is shown in Fig. 7. A penalty must be paid for the 
ease with which simulation has been achieved and that is that 
the states through have no physical significance. In 
order to determine some physical quantity (e.g., z, z', f, 6'), 
linear combinations of the states x^ through x^ are requ red. 
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COEFFICIENT VALUES FOR EQUATION (3.8). 
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Figure 7. Analog Model 

(Assumes no phase reversal through integrators) 

2 . Digital Model 

The digital model could easily be solved by choosing 
the states as being the same as those in the analog model, 
again, however, the chosen states have no physical meaning, 

A more desirable model is achieved by choosing the 
states as being heave, heave velocity, pitch and pitch velocity. 
These states can be obtained by algebraic manipulation of 
Eqs. (3.6) and (3.7). In their present form, each of these 
equations contains the transforms of the second derivative of 
both heave and pitch. This explicit interdependence, illustrated 
by the a^-b^ loop in Fig. 6, must be removed in order to c tain 
a workable model. This is accomplished by substituting Lc (3.6) 
into (3.7) and Eq. (3.7) into (3.6). 
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The states are chosen as indicated above, resulting 
in a state differential equation which describes the ship. 
This equation has the form 



x ( t ) = Ax ( t) +Bu ( t) +Cw ( t ) 



where the state vector is defined as 



x (t) 



the control vector is 



B ( t) 



and the wave-force vector is 



w(t) 



The coefficient matrices are 



A = 



0 

■0.2951 

0 

0.012 



■0 . 9362U 
0 



« e (t> 



S c (t > 



w 



M 



w 



■0.9362U +0.0207 
0 



-0.5742U -0.5742U -0.2961 



(3.9) 



0 

-1.2990 

1 

-2.0910 
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0 



0 



and 



B 



-0 . 2964U 2 
0 

-0.3706U 2 



C 



0 

0.4907 

0 

-0.0502 



-0.1624U 2 

0 



0.1432U 2 

0 

-0.0503 

0 

0.2865 



Equation (3.9) is readily simulated on a digital 
computer to give the ship's time response to control surface 
deflections and known wave forces and moments. 
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IV. CONTROL SYSTEMS FOR THE SEMISUBMERGED SHIP 



Optimal control theory [Ref. 5] was used to derive three 
controllers for the semisubmerged ship. The derivations of 
the control laws are based upon two assumptions: (1) the 

states in Eq. (3.9) are directly measurable, and (2) no 
measurement noise exists in the measurement of the states. 

A. PERFORMANCE MEASURE 

The performance measure gives a quantitative measure of 

the quality of response of the system under consideration. 

/\ 

With the exception of long waves (\>14) where the wave crest 
might impact with the upper platform and it would be more 
advantageous to follow the contour of the wave, it is desired 
to maintain the ship as stable as possible in the Earth-fixed 
coordinate system. This requires that the state vector be 
maintained close to the equilibrium point (origin) . The con- 
trol surfaces however, should not have excessive deflections 
which would run them into mechanical or electrical stops, nor 
should they overload the actuation system and result in loss 
of the control system. In view of the above, the requirement 
for the control system becomes: 

The control system must maintain the state 
vector, x(t) , as close to the origin as 
possible without an excessive expenditure 
of control effort. 

The performance measure then chosen is 
t f 

J = 1/2 J t [x T ( t) Qx ( t) +u T ( t) Ru ( t) ] dt (4.1) 

0 “ 



28 



where the final time t f is assumed to be fixed, Q is a real 
symmetric positive semi-definite weighting matrix and R is a 
real symmetric positive definite weighting matrix. It is 
assumed that the states and controls are not bounded and that 
the final state, x(t^), is free. No penalty has been placed 
on minimizing the deviation of the final state from its 
desired value, 0. It is desired to determine the control 
u*(t) that minimizes this performance measure. 

B. CONTROL LAW FOR THE SHIP IN A CALM SEA 

The system of equations defining the ship in a calm sea 
(w(t)=0) is 

x(t) = ^x+Bu(t) (4.2) 



The requirement to minimize the performance measure J, given 
by Eq. (4.1) leads to the well-known optimal control law 
[Ref. 5] 



u*(t) = -R -1 B T K(t)x(t) 

= F (t) x ( t) 



(4.3) 



where K(t) is the real symmetric positive definite solution to 
the Riccati equation 

K ( t) = -K(t)A-A T K(t)-Q+K(t)BR _1 B T K(t) (4.4) 



which satisfies the boundary condition 

K(t f ) = 0 

The optimal control law given in Eq. (4.3) is seen to be a time- 
varying feedback of the state vector x(t). The closed-loop 
system. state equations reduce to 

x ( t ) = [A+BF(t) ]x(t) (4.5) 
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when the optimal control law is used. Figure 8 shows this 
system in block diagram form. 



C. THE OPTIMAL CONTROL LAW FOR THE SHIP WHEN SUBJECTED TO 
WAVE FORCES 

The Hamiltonian for the forced system is 



(x(t) ,u(t) ,g(t) ,t) = 1/2x T ( t) Qx (t) +1/2 u T (t) Ru (t) 

+ 2 T (t)Ax(t)+£ T (t)Bu(t) (4 

+£ T ( t) Cw (t) 

where p(t) is a vector of Lagrange multipliers called costates. 
The theorem that the variation of the augmented functional 
fc f 

J = / { 1/2 [x T ( t) Qx (t) +u T (t) Ru ( t) ] 

a t 

0 (4.7) 

( t) [Ax (t) +Bu ( t) +Cw ( t) -x ( t) ] }dt 



equals zero on an extremal requires that the following neces- 
sary conditions hold for the state equation, costate equation, 
and control [Ref. 5]: 



x*(t) = Ax* ( t) +Bu* ( t) +Cw ( t) 



E*(t) = - 1^- = -Qx* (t)-A i p*(t) 



for all 
t= [t Q , t f ] 



0 = - |^ = Ru* (t)+B T p* (t) 



(4.8) 

(4.9) 
(4.10) 



where the asterisk denotes an extremal. R is a positive 
definite matrix, therefore Eq. (4.9) can be solved for the 
optimal control to give 

u*(t) = -R 1 B T £* (t) 4.11) 
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Figure 8. Linear Regulator Plant with Optimal Feedback Controller 



Substitution of Eq. (4.11) into Eq. (4.8) yields 



,-l„T 



x* ( t ) = Ax*(t)-BR B 2 *(t)+Cw(t) 

Equations (4.9) and (4.12) now form a system of 2n linear 
nonhomogeneous differential equations 



(4.12) 



x* (t) 
£* ( t) 



A 



-Q 



-BR B' 
-A T 



whose solution is 



X* (t f ) 




X* (t) 


— r 


= i(t f ,t) 








E* (t f } 




£* (t) 
— 



l 


X* (t) 


+ 


Cw ( t ) 




_£*(t)_ 




0 



(4.13) 



+ / $(t f ,T) 



Cw (t) 



di 



(4.14) 



where <£ is the state transition matrix of Eq. (4.13). If <j> 
is partitioned and the integral in Eq. (4.14) is replaced by 
a 2nxl vector 



/ <Ht f ,x) 



Cw (t) 



dx = 



f x (t) 



f 2 (t) 



(4.15) 



Eq. (4.14) can be rewritten as 

x*(t f ) = (t f , t) x* ( t) (tf • t) E* (t) + f]_ (t) 

£*(t f ) = ^21 (tf ' t) 2* ( t) + ^22 (tf ' t) E* (t) +f 2 (t) 

The boundary condition for the costates .is 



(4.16) 

(4.17) 



E (t f ) = 0 
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Therefore, Eq. (4.17) can be solved as 



E*(t) 


= - $ 2 ~ 1 (t f , t) $ 21 (t f , t)x* (t) 

(4.18) 




- ^2 2 1 (tf , t) £ 2 (t) 


Kalman (Ref. 6] has 


shown that the required inverse. 


^22 exists 

can be rewritten as 


for all t e [tg,tf] so that Eq. (4.18) 


B*(t) 


= K(t)x*(t)+s(t) (4.19) 


Equation (4.19) can 


now be substituted into Eq. (4.11) resulting 



in the optimal control law (Eq. (4.20)) similar to that in 



Ref. [7]. 




u* (t) 


-IT -IT 

= -R(t)B(t)K(t)x* (t)-R(t)B(t) s(t) 

(4.20) 




= F(t)x*(t)+v(t) 



Equation (4.20) indicates that the optimal control law contains 
a linear function of the system states, and also a command 
signal, v(t), that not only depends upon the system parameters, 
but also upon future values of the wave forcing function, w(t). 
This indicates that the optimal control law must be anticipatory 



in effect. 




The matrix K(t) 


and the vector s(t) can be obtained by 


differentiating Eq. 


(4.19) 



£*(t) = K(t)x*(t)+K(t)x* (t)+s(t) (4.21) 

and substituting Eqs. (4.8), (4.11) and (4.19) into Eq. (4.21). 
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g*(t) = K(t)x* (t)+K(t) [Ax*(t)+Bu*(t)+Cw(t) ]+s(t) 

= K(t)x*(t)+K(t) [Ax* (t) -BR _1 B T g* ( t) +Cw (t) ] +s (t) 

= K(t)x*(t)+K(t) [Ax*(t)-BR -1 B T (K(t)x* (t)+s(t) ) 



+Cw(t) ] +s (t) 



(4.22) 

However, substitution of Eq. (4.19) into the costate Eq. (4.9) 
also results in 



g*(t) = -Qx*(t)-A T (K(t)x*(t)+s(t) ) 



(4.23) 



Equations (4.22) and (4.23) can be equated which, after 
collecting terms, yields 



[K (t) +K (t) A-K (t) BR _1 B T K (t) +A T K (t) +Q]x* (t) 

' - *y —■» 

+[s(t)+A T s(t)-K(t)BR -1 B T s(t)+K(t)Cw(t) ] = 0 



(4.24) 



which must be satisfied for all x(t) and w(t) . This results 
in the matrix differential equation for K(t), 

K(t) = -K(t) A-A T K(t) -Q+K(t)BR _1 B T K(t) (4.4) 

which is the same as that given earlier for the calm sea 
condition, and an additional vector differential equation, 

s(t) = -[A T -K(t)BR _1 B T )s(t)-K(t)Cw(t) (4.25) 

The boundary conditions to be satisfied are 

K(t f ) = 0 

s(t f ) = 0 
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Equations (4.4) and (4.25) can now be integrated backward 
in time from t^ to tg and the resulting solutions can be 
stored and substituted into Eq. (4.20) to give the optimal 
control law 

u' (t) = -R _1 B T K (t) x* (t) -R _1 B T s (t) (4.20) 

= F (t) x* (t)+v(t) 

The controlled system, shown in Fig. 9, is now defined by the 
matrix differential equation 

x ( t ) = [A+BF(t) ]x(t)+Bv(t)+Cw(t) (4.26) 

D. A SUBOPTIMAL CONTROL LAW FOR THE SHIP WHEN SUBJECTED TO 

WAVE' FORCES 

The optimal control law given by Eq. (4.20) requires that 
the command signal v(t) be determined either by (1) integrating 
Eq. (4.25) backwards in time assuming the forcing signal w(t) 
is known, or (2) storing previously computed values of v(t) in 
a look-up table and selecting the one which corresponds to the 
appropriate time history of w(t). Either of the above mo-thods 
would require the use of either an existing onboard computer 
or a computer dedicated specifically to the control system for 
computation and storage. The suboptimal control system erived 
below was investigated to determine if this added complexity 
could be reduced or altogether avoided. 

Referring to the system equation given by Eq. (4.26), the 
objective is to have the state vector, x(t), approach 0. This 
requires that x(t) also approach 0. This can only occur if 
when x ( t ) - 0, 

Bv (t) = -Cw (t) 1 4 . 27 ) 
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In general, Eq. {A. 21) cannot be solved for v(t) since B~ 
does not exist; however here the first and third rows of 
both B and C are zero which allows two new matrices, B' and 
C', to be defined as 



B ' 



C' 



b 



b 

C 



21 b 22 

41 b 4 2 
21 C 22 



C 



41 




(4.28) 



It can be shown (see Appendix A) that 



B'v(t) = -C'w(t) (4.29) 

and since B' and C' are (2x2) matrices and B' ^ exists, Eq. 
(4.29) can be solved to yield 



v ( t ) = -B' 1 C'w(t) 

= Gw ( t) 



(4.30) 



The control law can now be written as 



u(t) = F (t) x(t) +Gw(t) 
and Eq. (3.9) becomes 

x ( t ) = [A+BF(t) ]x(t) + [BG+C]w(t) 



(4.31) 



(4.32) 



It can be easily shown that as defined, the matrix equality 



BG = -C 



(4.33) 
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holds (see Appendix B) so that if w(t) is exactly known, 

Eq. (4.32) reduces to Eq. (4.5), the system equation for calm 
seas . 

x(t) = [A+BF(t) ]x(t) (4.5) 

Therefore, if the suboptimal control law given by Eq. (4.31) 
is used, the computer capability required by the optimal 
control law is avoided. 

Reference [7] shows that the eigenvalues of the matrix 
[A+BF(t)] have negative real parts so that after initial 
transients, the system states will approach zero. Figure 10 
shows the block diagram for this system. 

E. CONSTANT STATE FEEDBACK 

The control laws given by Eqs. (4.3), (4.20) and (4.31) 

contain the time-varying feedback matrix K(t). This matrix 
was determined, as previously stated, by integrating Eq . (4.4) 

backwards in time. Although K(t) can be computed and stored 
prior to utilization of the above control laws (since the solu- 
tion of Eq. (4.4) does not depend upon the system states), 
the complexity of the control ' system would be increased over 
one which incorporated constant rather than time-varying 
feedback of the system states. 

Kalman [Ref. 6] has shown that if (1) the plant is completely 
controllable, (2) there is no cost penalty for final state vector 
deviations from the desired final state, and (3) A, B, R and Q 
are constant matrices, then K(t) approaches K (a constant 
matrix) as t^ approaches infinity. To determine whethe: the 
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Figure 10. Linear Regulator Plant with Suboptimal Feedback Controller for External 
Disturbances (w(t)) 



plant is completely controllable, the rank of the matrix 



E = [B AB A 2 B . . . A n 1 B] 



(4.34) 



must be evaluated. The plant is completely controllable if 
and only if E has a rank of n (in this particular case n is 
four). Equation (4.34) was evaluated and found to have a 
rank of four, so this system is controllable as anticipated. 

The ship meets all of the above requirements with the 
exception that t^ does not approach infinity (although it may 
be quite large). However, the solution to Eq. (4.4) using the 
NPS IBM-360/67 computer (Figs. 11 through 15) shows that the 
feedback gain matrix does reach steady-state constant values 
in a very short time ( — 3.1 sec. for a 2000-ton-displacement 
ship) and hence K(t) can be approximated by the constant matrix 
K to give the following control laws which are used in the 
remaining investigation: 



u(t) = Fx ( t) 
u*(t) = Fx ( t) + v ( t ) 

u(t) = Fx(t)+Gw(t) 



F can be obtained from 

-1 T 

F = -R B K 



(4.35) 

(4.36) 

(4.37) 



(4.38) 



where K is the steady-state solution to Eq. (4.4) and hence 
also satisfies the nonlinear algebraic equation 

0 = -KA-A T K-Q+KBR~ 1 B T K (4.39) 
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Figure 11. Time Histories of K 11 (t), K 1? (t) and 
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Figure 12. Time Histories of K 13 (t), K 23 (t), K 31 (t) and K 32 (t). 
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Figure 13. Time Histories of K 14 (t), K 24 (t) , K 41 ( t) and K 42 (t). Q - ±r 5 0.11. 
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Figure 15. Time Histories of K__(t) and 



V. SIMULATION PROCEDURE 



The NPS-360/67 computer was used to solve the ship's time 
response incorporating the control laws determined previously 
for various wave forcing functions. The computer program is 
given following the appendicies, and a flowchart is shown in 
Fig. 16. In order to solve Eqs. (4.25) and (4.32), w(t) and 
v(t) must be known. 

A. WAVE MODEL 

The wave forces were estimated using the properties of 
deep water Airey waves. The resulting equations are [Ref. 3] 

A A A 

z = Acos (-a ' t) +B sin (-a ' t) (5.1) 

w 

M = Ceos (-a ' t) +D sin (-a ' t) (5.2) 

w 

where 

A 

a' = non-dimensional wave' encounter frequency 

= ( 2 tt/A ) 1//2 - ( 2ttU/X ) cosi[> (5.3) 

= wave direction 

0° (following sea) 

= (5.4) 

180° (ahead sea) 

/\ 

and X is the non-dimensional wavelength of the ocean wave. The 
coefficients A, B, C and D are graphed versus X in Ref. [3]. A 
least squares polynomial curve fitting subroutine, LSQPL2 , was 
used to fit a sixth-order polynomial to each of the grap’- 3. 

/s 

Each polynomial is a function of X . The resulting polyn .nials 
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Figure 16. Computer Flowchart. 



are used in the computer program to simulate the wave forces 
and moments acting on the ship. 

B. COMMAND SIGNAL 

The optimal control law (Eq. (4.36)) requires the command 
signal 

v(t) = -R _1 B T s ( t) (5.5) 

where s(t) is the solution to Eq. (4.25). Equation (5.5) was 
solved for the three cases shown in Table VI. All three cases 
are for following seas (ij)=0°) which, as previously stated, is 
a sea condition which gives an undesirable ship response. The 
results for case three are shown in Figs. 17 and 18. It is 
seen that as expected for a sinusoidal forcing function, the 
command signal is also sinusoidal in form with the exception 



TABLE VI 

PARAMETERS USED IN SOLVING EQUATION (5.5) 



Case No. 


Q 


R 


A 

u 


/N 

X 




1 


I 


I 


1.0 


10.0 


0° 


2 


I 


0.51 


1.0 


10.0 


0° 


3 


I 


0.11 


1.0 


10.0 


0° 















I = Identity Matrix 



being where t approaches t^ and the boundary condition 

v (t f ) = 0 

must be met. 
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Since Eq. (4.25 must be solved forward in time and Eq. 

(5.5) can only be solved backward in time, it becomes necessary 
to have a model for v(t) to use in solving Eq. (4.25). The 
model that was chosen was of the form 



v 1 (t) = M^sin (-crx+y^) 
v 2 (t) = M 2 sin (-ct+y 2 ) 



(5.6) 



where the magnitudes, and M 2 , and phase angles, y^ and y 2 , 
were determined from the solutions to Eq. (5.5) and the final 
time deviations from sinusoidal form were ignored. The models 
for the three cases in Table VI are incorporated in the computer 
program as subroutine EST. 

An alternative approach for obtaining the command signal 
is to use a frequency approach instead of integrating Eq. (4.25) 
backward in time. If a new time variable, x , is defined as 



then 




dx = -dt 

and Eq. (4.25) can be rewritten as 

s ( T ) = [A T -KBR _1 B T ] s (t) +KCw(t) (5.7) 

which satisfies the initial condition 

s (0) = 0 

. —1 T 

As previously stated, [A-BR B K] is stable and since the 
transpose of a stable matrix is also stable, Eq. (5.7) is a 
stable equation in x (backwards in time t) . 
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The Laplace transform of s(t), after invoking the 
boundary condition, is given by 

S(s) = [sI-A T +KBR _ 1 B T ] 1 kcw(s) 



(5.8) 



If w(x) is a periodic function with a discrete Fourier spectrum 
W(f) , then the spectrum S(f) is given by 

S(f) = [j2TTfI-A T +KBR _1 B T ] 1 KCW(f) (5.9) 



and s_(t) can be obtained by taking the inverse Fourier trans- 
form of Eq. (5.9). This approach will hold for any general 
function w(t) since the initial condition of Eq. (5.7) is zero. 
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VI. RESULTS 

The three controllers derived in the previous chapter will 
be called (1) the "constant-gain controller" — given by Eq . 
(4.35), (2) the "optimal controller" — given by Eq. (4.36) 

and (3) the "suboptimal controller" — given by Eq. (4.37). 
Although the constant-gain controller was derived for the 
calm sea condition, its use was also investigated for the 
condition where w(t) / 0. Figures 19 and 20 show the ship's 
heave and pitch response to a following sea of wavelength 10 
when the ship's speed is 1.65 and no control system is used. 

A. CASES INVESTIGATED 

The following assumptions were made in the first five cases 
listed below. 

(1) The control surface deflections were 
instantaneous . 

(2) The wave forces and moments, w(t) , were non-zero 
and exactly known. 

(3) The state vector weighting matrix, Q, was equal 
to the identity matrix, I. 

(4) Contact of the underwater hulls and control 
surfaces with the air-water interface was ignored. 

(5) Initial conditions which were considered severe 

z(t Q ) = -0.5 

z(t Q ) = 0.0 
6 (t Q ) = 0.5 radians 
0(t Q ) = 0.0 radians/sec. 
were used (refer to Fig. 2) . 



53 





CM rH O rH CN 

i 1 



e 

<D 
p 
U ) 
>1 
cn 



o 

p 

P o 

C! rH 

0 

O ii 



0 


< r< 


£ 






o 


p 


O 


•H 




£ 


ii 


0) 




y) 




a 




0 


LO 


04 


KD 


U) 


• 


<d 


rH 




ii 


<D 




> 


<D 






Q) 




w 








<D 

u 

0 

tr> 

■H 

(*4 



54 



w 

£ 



u 




cn 



CN 



CN 

I 



55 



In the sixth case it was assumed that w(t) was not known 
exactly; the other assumptions indicated above remained in 
force . 

1 . Case I 

The reduction of the control weighting matrix R and the 
resultant effect on the ship's response and control surface 
deflections was investigated. 

2 . Case II 

The control surface deflections were constrained to 
+30° limits and the effect on the ship's ^response was determined. 

3 . Case III 

The ship's response using the Constant-Gain Controller 
was determined. 

4 . Case IV 

The ship's response using the Optimal Controller was 
determined . 

5 . Case V 

The ship's response using the Suboptimal Controller 
was determined. 

6 . Case VI 

The effect of imperfect estimates of w(t) on the ship's 
response was investigated using the Suboptimal Controller. The 
estimate of w(t) was perturbed by +15 percent of its true value. 

B. ' THE EFFECT OF VARYING CONTROL WEIGHTING MATRIX (CASE I) 

The values of the weighting matrix R in Eq. (4.1) determine 
how much emphasis is placed upon the expenditure of control 
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effort as compared to maintaining the state vector, weighted 
by Q, close to the origin. It is seen that as the values of 
R are made smaller compared to those of Q, more emphasis will 
be placed upon maintaining x(t) close to the origin. 

The ship's response in heave and pitch is shown in Figs. 
21 and 22 for a following sea with a wavelength of 10 and a 
ship speed of 1.0. The Optimal Controller was used with two 
different R matrices 



R 



1.0 0 
0 1.0 



and 



R 



0.5 0 

0 0.5 



The figures show, as expected, that heave and pitch are reduced 
as R is reduced. The corresponding increase in control effort 
is shown in Figs. 23 and 24 for the Elevator and Canard 
respectively . 



C. CONSTRAINED CONTROL SURFACE DEFLECTIONS (CASE II) 

A large canard deflection (see Fig. 24) was required to 
reduce the state vector to its steady-state value when R was 
reduced. The ship's response to the same conditions as in 
Case I above with R = 0.51 is shown in Figs. 25 and 26 where 
the control surface deflections were arbitrarily limited to 
+ 30°. The difference in state trajectories between the un- 
constrained and constrained control surface deflections s 
almost negligible; additional verification is provided the 
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Figure 23. Elevator Response Using Optimal Controller 
q = I, R = I and 0.51, U = 1.00, ^ = 0°, X = 10 
Control Surfaces Unconstrained 
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performance measure (see Table VII) . It was found that an 
additional reduction in the control weighting matrix to 



0.1 0 



R = 



0 0.1 



brought the state vector much closer to the origin but with 
unconstrained control deflections, produced unrealistic deflec- 
tion angles in the initial control period. The control surfaces 
were then constrained to the previous limits and the resultant 
differences in the state trajectories were negligible since 
the period of constraint was small compared to the total control 
interval [t^, t^] . As a result, all remaining cases have con- 
strained control surface deflections (+ 30°) and 



0.1 0 



R = 



0 0.1 



to bring the state vector closer to the origin. 



D. THE CONSTANT-GAIN CONTROLLER (CASE III) 

As previously stated, this controller was not derived for 
the case where w(t) ^ 0. However, since it is based upon feed- 
back of the system states, it will tend to reduce any disturbance 
of the state vector from its equilibrium point. Figures 27 
through 29 show the ship's response and control surface deflec- 
tions for a following sea with a wavelength of 10 when the 
ship's speed is 1.65. As expected, the state vector's devia- 
tion from equilibrium has been reduced (refer to Figs. 1' and 
20 for comparison) and the Constant-Gain Controller can e 
used without any knowledge of w(t). 
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Figure 29. Elevator and Canard Response Using Constant-Gam Controller 

q = I, R = 0.11, U = 1.65, ^ = 0°, X = 10 
Control Surfaces Constrained to + 30° 



E. THE OPTIMAL CONTROLLER (CASE IV) 

The ship's response using this controller is shown in 
Figs. 30 and 31. The ship's speed was 1.00 and the seas 
were following with a wavelength of 10. The control surface 
deflections (Fig. 32) are seen to be less than those of the 
Constant-Gain Controller as are the deviations of the state 
vector from equilibrium. This is expected since the Optimal 
Controller minimizes the performance measure J. 

F. THE SUBOPTIMAL CONTROLLER (CASE V) 

Figures 33 through 35 show the ship's response and control 
surface deflections for a following sea of wavelength 10 and a 
ship's speed of 1.65. The response and control surface deflec- 
tions for an ahead sea (^ = 180°) of wavelength 10 are shown 
in Figs. 36 through 38. It is seen that, as previously stated, 
the state approaches zero after the initial transients. 

Although the states are maintained at the origin after the 
initial transients, the control effort required is greater than 
that of the Optimal Controller and, as expected, the performance 
measure is greater (see Table VII) . 

G. THE EFFECT OF INACCURATE ESTIMATION OF WAVE FORCE (CASE VI) 
The response of the ship to imperfect estimates of w(t) was 

investigated for the case where only the estimates of the 
amplitudes were imperfect. It was assumed that the wavelength 
and phase were known exactly. 

The Suboptimal Controller was used to control the ship 
when the amplitude of w(t) was not exactly known. The r ip was 
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TABLE VII 



PERFORMANCE MEASURES 
t Q = 0 sec., t f = 80 sec. 



Controller 


Q 


R 


/\ 

u 


* 


A 


w (t) 


j 


None 


I 


0.11 


1.65 


0° 


10 


— 


9.328477xl0 4 


Optimal 

(unconstrained) 


I 


I 


1.00 


0° 


10 


Known 


8.5290623xl0 2 


Optimal 

(unconstrained) 


I 


0.51 


1.00 


0° 


10 


Known 


2 

6.5330823x10 


Optimal 

(constrained) 


I 


0.51 


1.00 


0° 


10 


Known 


6.5468393xl0 2 


Constant-Gain 

(constrained) 


I 


0.11 


1.65 


0° 


10 


Known 


3.5693990xl0 2 


Optimal 

(constrained) 


I 


0.11 


1.00 


0° 


10 


Known 


3.9204965xl0 2 


Suboptimal 

(constrained) 


I 


0.11 


1.65 


0° 


10 


Known 


3.2483721xl0 2 


Suboptimal 

(constrained) 


I 


0.11 


1.65 


180° 


10 


Known 


3.9344051xl0 2 


Suboptimal 

(constrained) 


I 


0.11 


1.65 


0° 


10 


Underestimated 
by 15 percent 


3.3556745xl0 2 



running in a following sea of wavelength 10 at a speed of 1.65 
The ship's response in heave and pitch are shown in Figs. 39 
and 40 where the amplitude of w(t) was overestimated and under 
estimated by 15 percent. It is seen that the estimate error 
has very little effect on the performance of the Suboptimal 
Controller . 
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VII. PRELIMINARY WORK ON WAVE ESTIMATORS 



The estimators stated below are possible candidates which 
can provide an estimate of w(t) required for the operation of 
both the Optimal and Suboptimal Controllers. 

A. AN INPUT-SIGNAL OBSERVER 

The estimator depicted in Fig. 41 consists of a model of 
the plant being run in parallel with the actual plant. The 
state equation of the model is given by 

e(t) = Ae (t) +Bu ( t) +G [x (t) -e (t) ] 

and the state equation of the plant is given by Eq. 
writing this equation, it has been assumed that all 
the system are measured. If the difference between 
and model states (x(t)-e(t)) is called a(t), then 

a(t) = x(t)-e(t) 

= A[x(t)-e(t)]-G[x(t)-e(t)] +Cw ( t) 

= [A-G] a (t) +Cw(t) 

The solution to Eq. (7.2) is 

alt) = e [*-S]t a( 0 )+/ t e [A-G)(t-T) Cw(T)dT 

0 

so that if the eigenvalues of G are chosen such that e ' - w 
dies out quickly as t increases, then c(t) can be closely 
approximated by 

t 

o ( t ) - / £ (t-x) Cw ( t) dx .7.4) 

0 



(7.1) 

(3.9). In 
states of 
the plant 



(7.2) 



(7.3) 
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where <£ is the state transition matrix of Eq. (7.2). Referring 
to Eqs. (5.1) and (5.2), it is seen that Eq. (7.4) can be 
rewritten as 



o (t) = / $(t-x) 

0 



a^coswT+a 2 sintoT 



$2COSooT+$2 s i nu ' T 



dx 



(7.5) 



The coefficients a^, a 2 , 3-^ and 3 2 are constants and can be 
taken outside the integral; thus, Eq. (7.5) can be rewritten in 
the form 

t t 

o.(t) = a n J <f> . 9 (t-x) costoxdx+a„J c|> . „ (t-x) sinioxdx 

l 1 0 1Z 1 0 

(7.6) 

t t 

+ 3-,/ <P ■ a (t-x ) cosooxdx+3-/ <j> . . (t-x) sinwxdx 
1 0 14 2 0 1 

(i=l, . . . , 4) 



The integrals in Eq. (7.6) can be precomputed for a time 
interval (0,t^) and stored in a matrix D(t^) so that Eq. (7.6) 
can be put in the form 

g(t 1 ) = D(t 1 ) 

If the state differences are now sampled and stored at times 
t^ , i=l, ..., k, then an estimate of the coefficients a^, a 2 , 
3^ and 3 2 can be obtained by determining the minimum-mean- 
square solution of 



a, 



a. 



3 



1 

32 



(7.7) 
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0 (t ± ) 




r ( )" 




V 




“i 


o(t 2 ) 

• 

• 

• 


= 


D (t 2 ) 

• 

• 

• 




a 2 

P 1 


= jD 


a 2 

S 1 


,2‘Vj 




e'V 

L, _ 




*2_ 




S 2. 



given by 



V 




*2 (tj^r 


a 2 

^1 


= J)* 


a (t 2 ) 

• 

• 

• 


_ e 2_ 




_2<t k ) 



where & # is the pseudo-inverse ofj). 



(7.8) 



(7.9) 



B. FAST-FOURIER TRANSFORM ESTIMATION 

An alternative approach is to Laplace transform Eq. (7.2) 
to give 



-1 

o (s) = s [sI-A+G] C (0) 

-1 

+ [sI-A+G] Cw ( s) 



(7.10) 



If w(t) is a periodic function, then the spectrum £(f) will be 



-1 



£(f) = [ j 2 it f I -A+G ] CW(f) 
and an estimate of W(f) can be determined by 
CW(f) = [j2nfI-A+G] l (f) 



(7.11) 



(7.12) 



where 2(f) can determined by taking the Fast-Fourier 
Transform (FFT) of the error between the predicted and actual 
responses. 
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c. STATE OBSERVER MODEL 

A third possibility is to expand the state equation given 
by Eq. (3.9) by defining four additional states. The matrix 
product 

0 



Cw ( t) = 



can be rewritten as 



Cw ( t ) = 



coswt+a 2 sinwt 



cosa)t+g 2 s i nwt 



0 0 0 

cosait sinoot 0 

0 0 0 



(7.13) 



0 

0 

0 



a. 



a. 



coswt sinwt 3. 



(7.14) 



The augmented state vector is now defined as 

z(t) 
z (t) 

0 ( t) 

x (t) = 6 ( t) 

“l 
a- 



(7.15) 



so that the state equation becomes 
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1 0 


0 


0 


0 










A 


1 

1 coswt 


sinwt 


0 


0 






B 


= 




! 0 


0 


0 


0 


x (t) 


+ 


— 






0 

1 

1 

1 


0 


cosut 


sinwt 






0 




0 




0 











(7.16) 



= A' (t)x(t)+B'u(t) 

The actual physical states (x^(t) through x^(t)) are measured. 

The observation matrix for the augmented system is 

H = [10] (7.17) 

If an observer, similar to Fig. 41 but with eight states and 
an observation matrix H, is now run in parallel with the 
actual plant so that the observer's state equation is 

e ( t) = A' (t)e(t)+B'u(t)+G(t) [He(t)-Hx(t) ] (7.18) 

then the difference between the observer and actual state 
equations can be written as 

c(t) = [A ' (t) +G ( t) H] c (t) (7.19) 

By choosing G(t) so that [A ' ( t) +G ( t) H] is stable, then the 
states of the observer, e(t), will approach the states x(t) 
and the error c(t) will approach zero. The states e 5 (t) 
through eg (t) will then be estimates of a^, a 2 , 8^ and 6 2 
respectively. 

Utilizing the above method, an example can be given for 
the first-order case where 

x(t) = -ax ( t) +acoso)t a>0 
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The states are augmented by defining x 2 so that 

x^(t) .= -ax^+(cos t)x 2 

x 2 (t) = 0 

or 

x(t) = A ( t) x ( t) 



The observer state equation is 

e (t) = A(t)e(t)+G(t)H[e(t)-x(t) ] 



where the observation matrix is 

H = [1 0] 



The matrix 

[A (t) +G (t) H] = 

must now be tested for stability which can be accomplished by 
defining a Lyapunov function [Ref. 8] 

V ( t) = 1/2 (a^ (t) +c 2 (t) ) (7.20) 

and its time derivative 

V (t) = a 1 (t)a 1 (t)+a 2 (t)a 2 (t) (7.21) 



-a+g- 



’ 2 ' 



COSCOt 



For this example, 

V(t) = (-a+g 1 ) c^ (t) + (coscot) a 1 (t) a 2 (t) +g 2 c 1 (t) a 2 (t) 

= (-a+g^) (t) + (coscot+g 2 ) a 1 (t) a 2 (t) 

where if the values of g^ and g 2 are chosen such that 



g x < a 
g 2 = -coscot 
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then Eq. (7.20) is positive definite and v(t)<0. The system 
is therefore stable and the estimator states will approach 
the plant states as t approaches infinity. 

If this approach is used, additional work will be required 
to ensure that G(t) can be selected so that [A ' ( t) +G (t) H] is 
stable for the ship. 
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VIII. CONCLUSIONS 



A. THE CONTROLLER 

The easiest controller to implement is the Constant-Gain 
Controller since, as previously stated, it in no way depends 
upon any knowledge (past, present or future) of the wave 
forces and moments. This controller would also be the least 
expensive since the other two controllers are essentially the 
Constant-Gain Controller with additional computed control 
signals. On the other hand, both the Optimal Controller and 
the Suboptimal Controller would be more expensive and complex 
than the Constant-Gain Controller but would maintain the state 
vector closer to the origin. 

The complexity of realizing the command signal depends on 
the nature of the wave forces and how they are estimated. If 
w(t) is sinusoidal in form and is estimated directly in the 
time domain, v(t) can be obtained directly for the Suboptimal 
Controller whereas the Optimal Controller requires either on- 
line or precomputed solutions of Eq. (4.25). In this case, 
the Suboptimal Controller is considered to be less complex and 
the "best" controller. If, however, w(t) is some general 
periodic function requiring the FFT estimator or if the fre- 
quency domain approach is used for the single sinusoidal w(t), 
then the difference in complexity between the Optimal and Sub- 
optimal Controllers is negligible and the choice of the "best" 
controller will depend upon how close it is desired to r ^intain 
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the state vector to the origin. In either case, the estimator 
output could be set to zero (if the estimates were bad) and 
both controllers would reduce to the Constant-Gain Controller, 
resulting in still satisfactory performance. 

B. RECOMMENDATIONS FOR FURTHER STUDY 

The assumptions used in this investigation concerning sea 
conditions, control surface deflections and measurements result 
in the recommendations given below. 

1 . Sea Conditions 

The constant height-to-wavelength ratio and single- 
frequency wave used in this investigation is not indicative 
of the true sea condition which is a spectrum of different 
height-to-wavelength ratios and frequencies [Ref. 2] . A more 
realistic model for sea conditions should be derived and the 
ship's response to this model investigated. 

2 . Control Dynamics 

The ship's response should be investigated for the 
case where the dynamics of the control surface actuation sys- 
tem are included. This investigation assumed an instantaneous 
control surface deflection capability. 

3 . Measurement Noise 

The ideal case of zero measurement noise was assumed 
throughout this investigation. Further investigation is recom- 
mended to include realistic measurement noise. 

4 . Wave Estimator 

A wave estimator should be investigated for incli sion 
in the ship model and the response determined using the stimator 
and measurement noises as noted above. 
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5 . Weighting Matrices 



Further investigation is recommended in the determina- 
tion of the values of the weighting matrices Q and R. The 
possibility of exceeding control surface or structural rates 
and the accompanying effects on personnel was not included 
in this investigation. If further investigation proves this 
to be a problem area, then new values for the state and control 
weights will have to be determined. 

6 . Wave Contouring 

The inclusion of a realistic sea condition model will 
change the point where upper-body contact with the wave occurs, 
previously stated as being wavelengths greater than 14. This 
point of contact should be determined and a controller that 
would allow the ship to follow the wave contour instead of 
driving through the wave should be utilized for this case. 

This would result in a linear tracking problem instead of the 
linear regulator approach used in this investigation. 
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APPENDIX A 



Here it is 



holds, where 



and from which 



The matrix 



B ' v ( t) 



and 



shown that the matrix equality 



B'v(t) = -C'w(t) 



(4.29) 



B' 



b 21 b 22 
b 41 b 42 



C' 



C 21 C 22 

5"41 C 4 2 _ 



the command signal, v(t), can be obtained as 



v (t) 



-B' 1 C'w(t) 



(4.30) 



product 



b 21 b 22 




v 1 (t) 


b 41 b 4 2 




v 2 (t) 




_ 



b 4 i v 1 (t)+b 42 v 2 (t) 
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C 21 C 22 




w x (t) 


C 41 C 42 




w 2 (t) 



C 21 W 1 ^ +C 22 W 2 ^ 

C 41 W 1 ^ +C 42 W 2 ^ 

Therefore, since it is given by Eq. {A. 21) that 

b 2 i v i (t) +b 22 v 2 = “ ( C 21 w i (t) +c 22 w 2 ( fc ) ) 



and 



b 41 v i<t)+b 42 v 2 (t) 
which can be written as 



b 21 V l (t)+b 22 V 2 (t) 



b 41 V l (t)+b 42 V 2 (t) 



or, 



B'v(t) 



“ ( C 4 1 w 1 (t) +C 42 W 2 ^ ^ 

C 21 W 1 ^ +C 22 W 2 ^ 
C 41 W 1 ^ +C 42 W 2 ^ 
- C'w(t) 
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APPENDIX B 



Here it is shown that the matrix equality 

BG = -C 

holds, which reduces the system equation to 

x = [A+BF (t) ]x(t) , 



(4.33) 



(4.5) 



a stable system whose state approaches zero as t approaches 
infinity, when the suboptimal controller is used. 

BG = B(-B ,- 1 C' ) 



21 



41 



0 

b 

0 

b 

0 

0 

0 

1 



22 



42 



= -BB' _1 C 


i 1 






b 42 


CN 

CN 

,Q 

1 




C 21 


C 22 


_b 41 


b 21 




i — 1 

a 

i 


C 42 



(b 21 b 42“ b 22 b 41 ) 



C 21 C 22 



C 4 1 C 42 



21 



41 



22 



42 



-C 
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SEMISUBMERGED SHIP PROGRAM 
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**IF N3 CONTROL SYSTEM IS USED, DATA CARDS II & 12 ARE NOT USED* 
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READ SHIP'S NGN-DIMENSIONAL SPEED 
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SUBROUTINE DTR AN ( X , Y »M , N) SH002980 

IS SUBROUTINE TAKES THE TRANSPOSE OF AN { M X N) MATRIX X AND RETURNS IT 
AN (N X M) MATRIX Y. (DOUBLE PRECISION) 
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DIFFERENTIAL "OJ'-T IONS FOP SOLVING THE CONSTANT K MATRIX. 

DERY ( 1) =-2.D0* ( A( 2» 1)* Y( 2)+A(4»l )*Y(4) )-Q(l»l)+P(2f2)*Y (2 )**2 SH00331 

* +(P(2»4)+P(4,2))*Y(2)*Y(4)+P(4,4)#Y(4)**2 SH00332 
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